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From the point of view of anisotropic flow this conference is quite remarkable. Elliptic 
flow has been fully recognized as an important observable providing information on the 
early stages of heavy ion collisions. An increased attention to anisotropic flow in recent 
years has resulted in a significant improvement in all three components of the field: tech- 
q ! niques and methods, analysis and presentation of the data, and theoretical understanding 
of the underlying physics. For a long time the anisotropic flow measurements suffered 
large systematic errors related to the so called "non-flow" contribution, the contribution 
due to the azimuthal correlations not related to the orientation of the reaction plane. 
At this Quark Matter conference, for the first time, we discuss the results mostly free of 
this uncertainty, using results obtained with 3-, 4-, and even more particle correlations. 
Elliptic flow up to p t ~ 10 — 12 GeV/c, the "wiggle" in directed flow, V2(pt) for a number 
of identified particles up to p t ~3-4 GeV/c, are only a few other results presented at this 
conference. 

The starting point for many recent developments was the Quark Matter '95 conference 
at Monterey. Then, the E877 Collaboration reported the first observation of anisotropic 
flow in ultra- relativistic nuclear collisions at the BNL AGS jlj], the discovery, which M. 
Gyulassy in his concluding remarks called the "major highlight of the year" . That analyses 
was done with the use of a new approach ||, the very notations of which, like v\ and 
t>2, are a common language now. The E877 Collaboration also presented the results from 
the first attempts to perform HBT interferometry of the anisotropic source, one of the 
important directions in today's HBT analysis. At the same conference, J.-Y. Ollitrault, 
foreseeing future precise measurements, emphasized the importance of the understanding 
and measuring of the non-flow correlations. It was pointed out that flow is a collective 
phenomena and in order to disentangle it from other contributions to particle azimuthal 
correlations one has to exploit multi-particle correlations (cumulants). 

Anisotropic flow is defined as azimuthal asymmetry in particle distribution with respect 
to the reaction plane (the plane spanned by the beam direction and the impact param- 
eter). It is called flow for it is a collective phenomena, but it does not necessarily imply 
hydrodynamic flow. It is commonly characterized using Fourier decomposition of the az- 
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imuthal distributions ||. Then the first harmonic Fourier coefficient, v\, describes directed 
flow, and the second harmonic coefficient, v 2, corresponds to elliptic flow; non-zero higher 
harmonics can be also present in the distribution. 

The two reasons for anisotropic flow are the original asymmetry in the configura- 
tion space (non-central collisions !) and rescatterings. In a case of elliptic flow the 
initial "ellipticity" of the overlap zone is usually characterized by the quantity e = 
((y 2 — x 2 ))/((y 2 + x 2 ), assuming the reaction plane being the xz-plane. With the system 
expansion the spatial anisotropy decreases. This is the reason for the high sensitivity of 
elliptic flow to the evolution of the system at very early times f|||, 2-5 fm/c, of the order 
of the size of the system. 



2. Directed flow 

Three major results reported at this conference regarding directed flow were: the un- 
derstanding of the role of momentum conservation ||[7| (due to lack of space I do not 
discuss it here, for details, see ||), the use of 3-particle correlations in order to suppress 
the non-flow contributions 0, and the first results on directed flow of protons close to 
mid-rapidity ||. The latter is related to earlier predictions of the so called "wiggle" in 
the dependence of vi(y) p|,^0|. 



Three-particle correlations, discussed in the talk of Borghini et al. 0, relies on the ob- 
served strong elliptic flow at mid-rapidity and permits directed flow measurements, which 
are much less affected by non-flow correlations. If in the usual two-particle correlator 
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vl + S 1 , (1) 



the non-flow contribution, S\ ~ l/N to t, where N to t is the total number of particles, then 
the 3-particle correlation approach gives 

( e ^i+<fe-2<« w V 2 V2 + 5 2_ (2) 

The relative suppression of the non-flow contribution in Eq. 2 compared to Eq. 1 is of the 
order of l/(v2N to t), which for large elliptic flow could be significant. Mathematically, one 
can trace that in Eq. |2], the contribution of any correlations projected onto the reaction 
plane and the plane perpendicular to the reaction plane enter with opposite signs. The 
non-flow contributions, by definition, are the correlations not dependent on the relative 
orientation with respect to the reaction plane, and as such just cancel out from the final 
result. 

The so called "wiggle" in directed flow, the three times change in the sign of directed 
flow in the midrapidity region, had been predicted a few years ago. The origin of this 
wiggle is somewhat different in different models, but any of the explanations would be very 
intriguing and interesting. In the "third flow component" picture ||, it results from the 
expansion of the initially tilted disk of the compressed matter. It appeared that the effect 
is present only if QGP is produced. Thus the wiggle might be a signature of the QGP 
production. The second explanation of the wiggle relies on the non-complete stopping 
of nucleons during the collisions, which results in the correlation of the position of the 
nucleon in the transverse plane and their rapidity. Combining this with a subsequent 
radial expansion one gets the wiggle fl(|. The first results have been presented by the 
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Figure 1. Directed flow of pions and protons at beam energies of 4CL4 and 160^4 GeV. 
Note the different scales in the two proton flow graphs. 



NA49 Collaboration || (see Fig. 1): One can see a clear indication of the change of sign 
of directed flow of protons as a function of rapidity for the most peripheral events. It 
would be extremely interesting to see the analogous results from experiments at RHIC. 
Detailed model predictions for the energy dependence of directed flow of both pions and 
baryons would also be very interesting. 

3. Elliptic flow 

There have been many new results and ideas presented at this conference. I will only 
briefly mention the most important (from my point of view). 

3.1. Comparing the results from different experiments 

Anisotropic flow results at this conference have been presented by many collabora- 
tions [P,[7|, pT1 , P^| , P^ ] . As always in such a case the question appears if these results are 
consistent. There are obvious difficulties to perform a fair comparison: different collab- 
orations use different centralities. Nevertheless, I would conclude that the data where 
comparable are consistent (for the figures see the slides of my talk). 

The second question often asked during the conference concerns transverse momentum 
dependence of elliptic flow at different collisions energies. The integrated elliptic flow is 
higher at RHIC than at SPS. The question is if v<i as a function of pt become steeper or if 
the observed increase in integrated flow is mostly due to the increase in mean transverse 
momentum? To answer this very interesting question is more difficult, as the slope of 
v z(Pt) strongly depends on centrality and one does have to have the same centralities in 
different experiments for comparison. With the existing data, however, I would conclude 
that in the central collisions the slope of i^G *) at SPS energies looks significantly smaller 
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than at RHIC; in mid-central collisions the difference is not that large. (I discuss the 
centrality dependence of v-i in more detail below.) 

3.2. Flow and non-flow from multi-particle correlations 

Anisotropic flow is a multi-particle phenomena. It means that if one considers many- 
particle correlations instead of two-particle correlations, the relative contribution of non- 
flow effects (due to few particle clusters) would decrease. Considering many-particle 
correlations, one has to subtract the contribution from correlations of the lower-order 



multiplets and use cumulants instead of simple correlators [[33]. For example, correlating 
four particles, one gets (itj = e*^*) 

(?Vl*V2<, 3 <,4> = + 2 ■ 2 ■ V 2 J n + 28% , (3) 

where 5 n stands for the non-flow contribution to two particle correlator: (M n l u* 2 ) = 
v 2 +5 n . Subtracting from the expression (§) twice the square of the two particle correlator, 
one is left with only the flow contributions 

where the notation ((...)) is used to denote the cumulant. A very elegant way of calcu- 
lating cumulants in flow analysis with the help of the generating function was proposed 
in p3[ . The NA49 Collaboration presented elliptic flow results from 4-, 6- and even 
8-particle cumulants [0]. Remarkably, the results from all higher order cumulants agree 
well, indicating that even 4-particle cumulants already suppress the non-flow contribution 
to a negligible level. The use of the multi-particle correlations permits one to measure 
the non-flow contribution and check its multiplicity dependence. The latter is expected 
to be roughly inversely proportional to the multiplicity, thus reflecting the small cluster 
origin of such correlations. Indeed, the result reported by the NA49 Collaboration are in 
agreement with a 1/N dependence. 

The STAR Collaboration presented results from 4-particle correlations at ^/snn =130 GeV 



and ^snn =200 GeV [I7[], Fig. 2. Transverse momentum dependence of the non-flow (rel- 
ative) contribution was also discussed |15[] . The results are consistent with weak increase 
in the relative non-flow signal from 15% at low p t to about 20% at p t ~ 4 GeV/c (at 
v /J/v r /v =130 GeV). Finally, note a possible increase of the relative non-flow contribution 
with the collision energy, from 7-10% at SPS |?[] to about 15% at ^/snn =130 GeV and 



~20% at tJsnn = 200 GeV at RHIC Both effects are consistent with the increase 
in the contribution from hard processes. 

3.2.1. Flow fluctuations 

Very good statistics results reported at this conference become sensitive to another 
effect usually neglected in flow analysis, namely, event-by-event flow fluctuations. The 
latter can can be due to two different reasons: "real" flow fluctuations - fluctuations 
at fixed impact parameter and fixed multiplicity (see, for example, [[HJ), and impact 
parameter variations among events from the same centrality bin in a case where flow does 
not fluctuate at fixed impact parameter. Note that these fluctuations affect any kind 
of analysis, including the "standard" one based on pair correlations. The reason is that 
any flow measurement is based on correlations between particles, which are sensitive only 
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to certain moments of the distribution in V2- In the pair correlation approach with the 
reaction plane determined from the second harmonic, the correlations are proportional to 
v 2 . Averaging over events gives (v 2 ), which in general is not equal to (v ) 2 . The 4-particle 
cumulant method involves the difference between 4-particle correlations and (twice) the 
square of the 2-particle correlations. It is usually assumed that this difference comes from 
non-flow correlations. Note, however, that this difference ((f 4 ) — (v 2 ) 2 7^ 0) could be due 
to flow fluctuations. Let us consider an example where the distribution in v is flat from 
v — to v — f max and there is no non-flow contribution. Then, a simple calculation would 
lead to the ratio of the flow values from the standard 2-particle correlation method and 
4-particle cumulants as large as (v 2 ) 1 ^ 2 / (2(v 2 ) 2 — (i> 4 )) 1//4 = 5 1//4 ~ 1.5. Thus comparing 
any theoretical model with the data one should take into account flow fluctuations if they 
exist in the model. 
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Figure 2. STAR Preliminary results |I7 on Figure 3. PHOBOS preliminary results fl~3|| 
centrality dependence of elliptic flow from on pseudorapidity dependence of elliptic 
4-particle cumulant analysis. flow at ^s^n =200 GeV. 



3.3. Pseudorapidity dependence 

The PHOBOS Collaboration is still the only one which reports elliptic flow in a wide 
pseudorapidity range, Fig. ^. This pseudorapidity dependence is still a challenge for 
the hydrodynamic models [|ll| , which predicts rather weak pseudorapidity dependence of 
elliptic flow. As was pointed out at the conference such behavior of v 2 in the hydrodynamic 
model could be related to the boost invariant initial conditions used in the calculations. 
Note also that PHOBOS results for elliptic flow at ^/snn =130 GeV and 200 GeV are 
very similar, which is in some disagreement with the STAR results. 

3.4. Elliptic flow at large p t 

An obvious interest in elliptic flow at large transverse momenta is due to a possible 
strong signal as a results of high p t parton energy loss [19,20]. Shuryak in his recent paper 



2~Tf argued that the signal observed by the STAR Collaboration (only elliptic flow results 



from 2-particle correlations were known at that moment) could be probably too strong. 
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He has calculated elliptic flow for the limiting case of very strong absorption resulting in 
predominant surface emission. Using his idea I arrived for the maximum flow values as a 
function of impact parameter 

V2,max = sin(2a)/(6a) , where a = arccos(6/(2i?)). (5) 



Comparing this formula to the STAR data [TJ] from 4-particle cumulants one finds rea- 
sonable agreement. 



At this conference, the STAR Collaboration presented the results |TJj on v 2 (pt) up to 



p t ~ 12 GeV/c (from 2-particle correlations). It shows that elliptic flow is large at least 
to pt ~ 6 — 8 GeV/c. Elliptic flow is maximum at pt ~ 3 GeV/c and probably decreases 
very slowly up to the highest momentum measured. To discuss reliably any trends in 
the region of p t > 6 GeV/c, however, requires multi-particle correlation measurements. 
Model calculations are also needed. 

Another interesting fact concerning flow at high p t has been reported by the CERES 



Collaboration [TJJ. According to their data, the saturation of v 2 (pt) at SPS energies is 
quite similar to the one at RHIC, although to make a definite conclusion one needs more 
carefully comparison. 

3.5. Baryon and meson elliptic flow 

At this conference both the STAR and PHENIX Collaborations presented very interest- 
ing results on p t dependence of elliptic flow of baryons and mesons. Although the results 
from both experiments have rather large error-bars, both agree that baryon elliptic flow 
saturates at higher p t and larger v 2 values than meson elliptic flow. From my point of view 
it could be due to the (partial) particle production via (constituent) quark coalescence at 
moderate transverse momentum. Then, in this region: 

^«[|^(P,«W2)] a , and ^ oc fe(p, - p B /3)]\ (6) 
d 6 p M d 6 p q d A p B d A p q 

resulting in larger saturation values of v 2 for baryons as well as the saturation point being 
at higher transverse momentum. 

Note that Eq. ^| is valid only in a rather limited region of transverse momentum. At 
very low momentum probably most hadrons are formed via coalescence. Then their 
elliptic flow values should be similar to the elliptic flow of constituent quarks (I do not 
consider here the possible dependence on the mass of the particles). At really high p t , 
however, hadrons are formed mostly due to parton fragmentation. There, the elliptic 
flow of hadrons should be determined by the elliptic flow of the fragmenting partons and 
there should be no significant difference in baryon and meson flow (due to this mechanism 
and neglecting the difference in the parton fragmentation functions into baryons/mesons). 
Note that the coalescence effect, if present, would also lead to an increase of the relative 
baryon to meson ratio in the same pt range. 

The centrality dependence of the ratio v 2t b ar yon(pt) / V2,meson(pt) would be extremely in- 
teresting. If the effect of baryon flow being larger than meson flow is indeed due to 
coalescence, one would expect that at higher centralities the crossings would appear at 
somewhat higher transverse momentum. A detailed comparison of the transverse momen- 
tum and centrality dependencies of elliptic flow of baryons/mesons and their yields could 
also help to understand the nature of both effects. 
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Figure 4. Upper panels: PHENIX Preliminary data on elliptic flow of protons/anti- 
protons and pions. Lower panels: STAR preliminary data on elliptic flow of A's and K° 
and protons and pions. 



3.6. Centrality dependence 

In the hydrodynamic limit elliptic flow is approximately proportional to the original 
spatial ellipticity of the nuclear overlapping region [^,^J2^], t> 2 oc e. In the opposite 
limit, usually called the low density limit [p^j^6[ , elliptic flow depends also on the particle 
density in the transverse plane: v 2 oc e dN/dy /S, where S is the area of the overlapping 
zone. The comparison of the results on elliptic flow from the point of view of particle 
density in the transverse plane was first done in ||26|| . In this approach, the transition to 
deconfinement would lead to some wiggles in v 2 /e dependence, ("kinks") [p^, p5| , p^] ). 

Fig. || shows the recent (^/snn = 130 GeV) STAR results on elliptic flow from 4-particle 
cumulants |I4] |, preliminary NA49 (also from 4-particle cumulants), and E877 results. 
New measurements at full RHIC energies of ^/snn = 200 GeV reported by STAR |17j are 



also shown. For this Figure I have rescaled down these new data points by a factor of 1.06 
to account for the change in mean transverse momentum (in 200 GeV data set the low 
transverse momentum cutoff is 150 MeV/c compared to 75 MeV in 130 GeV data). Note 
also that SPS data do not have a low p t cutoff. (There are possible systematic errors in 
Fig. 8 of the order of 10-20% due to uncertainties in centrality measurements. However, 
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these uncertainties do not alter the general trend). Fig. [| shows that at high particle 
densities the data reach hydrodynamical limits. With relatively large error-bars for the 
most central collisions it is difficult to draw any conclusion about the saturation at the 
hydrodynamical limit. As for the hydrodynamic limits, they are also somewhat model 
dependent. In some cases (depending on the latent heat and equation of state) in the 
model one can observe an almost linear dependence of V2 on particle density. 

A closer examination of the plot in the region of 7-15 particles/fm 2 shows that there 
is no scaling in this region. (This is probably related to the use of e calculated using the 
specific weight according to the number of wounded nucleons. One should also take into 
account slightly different transverse momentum regions used in different analyses.) Note 
that the data from both RHIC and SPS results show a rather flat region with some rise 
at density 15 fm~ 2 . One could speculate that such a behavior could indicate a change 
in the physics of rescatterings. Noteworthy that the color percolation point discussed by 
Satz [^HJ is just in the same region (shown by an arrow in figure). Very interesting that 
the WA98 Collaboration analyzing the high transverse momentum photon flow also 
observes a non-smooth behavior at the same place. 

Another indication of non smooth flow dependence on the particle density in this region 
can be seen in the figure showing an elliptic flow excitation function || . In this figure both 
NA49 and NA45 data points show very little increase in the average flow from 40A GeV 
to 160A GeV data, while there is a significant increase between highest AGS point and 
lowest SPS, and between SPS and RHIC.Q 
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3.7. Models 

It has been mentioned many times that at RHIC energies the hydrodynamic description 
of the results is quite successful. However, a few problems, like the description of the 

1 In the poster of H. Sako, CERES Collaboration, one can also notice an interesting increase in multiplicity 
as well as mean transverse momentum fluctuations at beam energy of 40^4 GeV. 
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pseudorapidity dependence [18] still remain. The hydro dynamic model has been discussed 
in a separate talk [31| at this conference and many details can be found there. 

At this conference there also was presented a model [32] with totally hadronic rescat- 
terings starting at zero time, which also describes many features of the data. 

AMPT (A Multi Phase Transport) model [30| also describes the data very well in the 
so-called "string melting" scenario with quark-quark cross sections of about 5 mb. More 
detailed examination of the "string melting" scenario reveals that it is nothing else than 
what one would do in a case of a constituent quark picture. Namely, the number of 
quarks in the system is calculated as the number of quarks in the final hadrons. The 
final formation time (after which a quark can interact) is also taken into account. The 
cross section of 5 mb also does not look unusual - it would be close to the total NN cross 
section divided by 9. Note that the AMPT model also quite successfully describes the 
HBT radii. This model definitely deserves further study. 



4. Conclusions 

I will summarize the most important results just by enumeration: 

- the "wiggle" in proton directed flow has been observed at SPS 0. 

- non-flow contribution is better understood and measured both at SPS and RHIC J7|,[l7 

- at RHIC V2(pt) saturates at values close the ones in the surface emission picture. 



- v 2 of baryons exceeds v 2 of mesons fromp t about 2 GeV/c till at least 3-4 GeV/c [12| . |T5 | . 

- At RHIC elliptic flow persists up to p t ~ 8 — 12 GeV/c 

- the saturation of elliptic flow at SPS probably not very different than at RHIC JTT[] 



- The dependence of elliptic flow on pseudorapidity is rather strong. It drops about 2 
times at about rj = 3 ][H|. (Still a challenge for the hydrodynamic description.) 

- v-ije is roughly proportional to particle density in the transverse plane (dN/ dy) / S reach- 
ing the hydrodynamic model predictions in central collisions. 



In this summary I also want to mention a few important measurements for the future: 

- we have to measure non-flow correlations in pp collisions and compare them to the ones 
measured in nuclear collisions. 

- we need to check for non-flow contribution at high transverse momenta. 

- centrality dependence of baryon/meson elliptic flow would be very interesting. 

- the collision of lighter system (e.g. Cu+Cu) and (multi) strange particle flow will be very 
important to understand the origin of the observed collectivity. 

- different RHIC collaborations have data taken at ^/snn =20 GeV. It would be really 
interesting to analyze these data and compare with SPS results. 

- directed flow at RHIC will be probably the next result from RHIC... 

- the percolation point region should be looked at in great detail. 

- finally, we need a detailed comparison (at a level of a few percent) of the results from 
different experiments; for that we need the results corresponding to the same centralities. 

Acknowledgments. I thank the organizers for the invitation, all the speakers for their 
help providing me their data, and many of my colleagues, and in particular Art Poskanzer, 
for numerous very helpful discussions. 
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